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ATTAINABLE CIRCULATION OF AIRFOILS IN CASCADE 
By Arthur W. Goldstein and Artur Mager 


SUMMARY 

From consideration of available information on boundary- layer 
behaviour, a relation among profile thickness, maximum surface 
velocity, Reynolds number, velocity diagram, and solidity is 
established for a cascade of airfoils immersed in a two- dimensional 
incompressible fluid flow. Several cascades are computed to show 
the effect of various cascade design parameters on minimum required 
cascade solidity. Comparisons with experimentally determined blade 
performance show that the derived blade loadings are equal or higher 
for moderate flow deceleration and somewhat lower for large decel- 
eration. Blades with completely laminar flow appear practical for 
impulse or reaction blading. 


INTRODUCTION 

In the design of compressor and turbine blading, the choice of 
blade sections and spacing is usually made on the basis of experi- 
mental results from several cascades. An enormous amount of data 
is required to cover the entire range of inlet and outlet angles, 
Reynolds numbers, blade sections, spacings, and mavi mum surface 
velocities. 

The relations among the blade circulation and its thickness, 
Reynolds number, maximum surface velocity, and spacing were deter- 
mined at the NACA Lewis laboratory in a qualitative fashion for 
two-dimensional, incompressible flow. This problem has been studied 
empirically (references 1 and 2) without any attempt to relate the 
problem to the basic determining factor - the boundary layer. The 
analysis presented herein establishes such a relation on a more 
rational basis. 

As the gas flows through a cascade of airfoils, on the suction 
surface of a blade the inlet velocity is raised to some mAvimim 
local value and then decreased until, near the trailing edge, its 
value is in the neighborhood of the exit velocity. A blade cascade 
will operate in a compressible fluid with less likelihood of local 
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shocks, resultant, losses, and choking if the design is such that 
the may 1 mum surface velocity is maintained as lov as possible. The 
preassigned, lov maximum surface velocity should extend over as 
large a portion of the suction surface as possible before diffusion 
in order to raise the circulation to a maximum. Because the- veloc- 
ity must decrease near the trailing edge froim its maximum value to 
a value close to the exit velocity, the maximum possible rate of 
diffusion will, in effect, determine the arc length on which the 
suction- surface velocity is kept at its maximum. The blade cir- 
culation attainable is then determined by the maximum diffusion 
rate possible. 

In dealing vith incompressible flow, the same general type of 
velocity distribution is specified herein, and the diffusion rate 
is determined by several conditions to avoid separation of flow 
from the blade surface. Once the possible velocity distribution 
on the suction surface is determined, an approximate computation 
is then used vith a specified velocity diagram and blade thickness 
to compute the pressure-surface-velocity distribution, blade cir- 
culation, and spacing. 

The results reported herein are not to be interpreted as giv- 
ing an exact velocity distribution of airfoil sections for appli- 
cation. In the first place, the assumption of no separation in 
some cases reduces obtainable lift and increases the blade solidity 
and the surface for friction, which increases skin friction loss. 

If a small region of separated flow is permitted, the increase in 
form drag may be insufficient to increase total drag in view of the 
increase in blade lift and decrease in blade number and skin fric- 
tion. Another limitation is possible in applicability of the 
boundary- layer equations used in prescribing the suction-surface- 
velocity distribution. The velocities deduced from these equations 
are not conventional in that the form factor of the boundary layer 
(ratio of displacement to momentum thickness) rapidly rises to the 
maximum permissible value and then remains nearly constant; whereas 
with conventional velocity distributions, it slowly increases and- 
then very rapidly rises near the separation region. Thus the 
empirical equations describing the development of the boundary- 
layer form factor may not be accurately applicable for such velocity 
distributions. In order to demonstrate the several solutions that 
might be applicable, the surface velocities were prescribed accord- 
ing to three rules. Gamer's equations (reference 3) were used 
with the form factor first rising rapidly and then remaining nearly 
constant with limiting values of either 1.85 or 2.14. Also a dif- 
fusion rate was chosen equal to 0.885 of the value prescribed by 
Kalikhman (reference 4) for separation. 



NACA TO 1941 


3 


A final consideration lies in the approximate nature of the 
procedure for computing the lower surf ace- velocity distribution. 

In this process, the effect of thickness is assumed to result in a 
surf ace- velocity component, which varies as a second-degree para- 
bolic function of the surface arc length. Consequently, the 
results apply to a special family of thickness distributions, which 
may not be the optimum for the prescribed suction- svirf ace- velocity 
distribution. Moreover, the estimate of the effect of thickness 
is inexact so that the estimate of the airfoil-thickness distri- 
bution will differ somewhat ' from the actual thickness of an airfoil 
with the prescribed surface velocities. 

Although the absolute values of blade circulations are inexact, 
some very definite trends are nevertheless demonstrated, particu- 
larly with respect to maximum surface velocity, Reynolds number, 
and blade thickness. 


SELECTION OF VELOCITY DISTRIBUTION 
ON SUCTION SURFACE OF BLADE 
General Considerations 

As mentioned in the Introduction, the velocity distribution on 
the suction surface to be investigated will consist of a velocity 
that is supposed to rise to it6 maximum value in a negligibly short 
distance from the leading- edge stagnation point and to maintain that 
value as far back along the airfoil as possible. The velocity must 
then decrease from the maximum value to the value near the trailing 
edge that differs slightly from the downstream velocity. A decel- 
eration with a maximum safe velocity gradient without flow separa- 
tion between the maximum velocity and a final velocity is required. 
In order to find the suction- surface -velocity distribution, it is 
necessary to investigate the development of the boundary layer, 
which determines the allowable rate of diffusion. The general 
shape of the velocity-distribution curves is shown in figure 1, 
which also Indicates some of the nomenclature. (All symbols used 
are defined in appendix A.) 


Laminar Boundary Layer 

Region of constant velocity. - For the region of constant 
velocity U(x) = U^^ , the equation for the momentum thickness of 

the laminar boundary layer is 
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d 


0.664 



( 1 ) 


where 

6 momentum thickness of boundary layer 

v kinematic viscosity of gas 

z arc length on airfoil suction surface 

U suction- surface velocity Just outside boundary layer - 

The subscript maz Indicates maximum value of surface velocity. 
The Reynolds number of the boundary layer based on momentum thick- 
ness is defined as 


R = U0/v 


( 2 ) 


The form factor H is defined as 

H £ &*/6 (3) 

where 5* is the displacement thickness of the boundary layer. In 
the region of constant velocity, H = 2.614. 

Region of diffusion. - For diffusion with a laminar boundary 
layer, the approzimate method of Loitsianskii (reference 5) was 
used. According to that system of equations, separation will occur 
if the proportional change in velocity dU/U per unit Of arc 
length dz equal to the generalized momentum thickness L reaches 
a value of -0.08884. That is 




L dU 
U dz 



where 


L £ 6R = e 2 ^ ( 5 ) 

must be greater than -0.08884. A value of = -0.06618 was 
arbitrarily selected to give a high diffusion rate without intro- 
ducing laminar separation. Loitsianskii’ s equation for boundary- 
layer growth then becomes 
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°L (£_ U\ 
i du _ F V a x 2 / 

(S) 


where F, a function of Is a constant equal to 0.8517, for 

\ - 0.06618. 

This equation is integrated to give 


dU 

dx 


^ = KU 


-F/Dl 


By substituting for du/dx from the definition of the con- 

stant of integration K is evaluated to give 



(6a) 


where the subscript 7 indicates the value at the beginning of the 
laminar diffusion process. By differentiating equation (6a) and 
eliminating dU/U by means of equation (4), there Is obtained 


dx dU °L . dL 
^ L = U = Dl+F L 

which is Integrated to 

i » 1 + - 1 * 0.7855 

Equations (l) and (5) are used to relate the constants 
and Iq 


(6b) 


*1 


From equation (5) 


I q = 0.441 Xj 


2 
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or 


6_ 

9, 


ii' 

A U , 


1 

,2 




F 

2(F+D l ) 


'_L_ 

A. 


,0.5421 


also 


(6c) 


E_ 

*1 


F+21^ 

^(F+Dl) 


B, "I I, 


_L_ 

A 


,0.4579 


(6d) 


The velocity potential cp, defined by 


cpS Udx 

do 


is found from 


cp - cp 


h u z h u z 


Udx = 


F+21^ 


L 

h 


F+2Dr 


vH - II 


= 1.390 


'_L_ 

A. 


,0.9158 


-1 


(6e) 


The constant cp^/u^ Lj is computed from the equations for the non- 
diffusing section giving 

, 

2.268 (6f ) 

♦ 

The form factor of the boundary layer H in the region Dr ■ -0.06618 
has a value of 3.214 (reference 5). 

In applications, it is convenient to know the position Xj in 

percent of the total length of the suction surface Z. From equa- 
tion (6b), 
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h. = h * A. (h ' l ) (7) 

The subscript f applies to the final value at the trailing edge 

x i l 

of the blade. In terms of velocities, with ■— = ■■■ ■ , 

lo vJ *441 


Then 


X_ 


2 - 268 (A) 



i 



u 


*1 

1 


X “ 


F+Dl 

i + 0 -i 41 

fk 

)Vx 

F+D L 

[W, 

) 


(7a) 


( 8 ) 


In figure 2, the variation of U/Uj and. r/Ej with x-Xj/Lj 

is shown for diffusion with the laminar boundary layer. In apply- 
ing these curves for any airfoil, the scale of the diffusing region 
must be adjusted for the value of Lj that is equal to 0.441 Xj. 

In figure 3, the curve for the entire suction- surface potential is 
shown with X as the unit of length rather than Lj . 


Transition 

As the air flows over the blade surface, the thickness of the 
laminar boundary layer will Increase until a point is reached where 
the character of the boundary- layer flow will change from 1 ami war 
to turbulent. The boundary layer will then develop according -to 
different laws and will therefore also sustain a different pressure 
gradient without separation. 

Transition occurs at a certain value of R^, which will depend 
on the turbulence intensity and scale, the pressure gradient, the 
roughness and the curvature of the surface, and on heat transfer 
through the boundary layer. An value of 250 has been selected 
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for this analysis in consideration of the high level of free-stream 
turbulence in turbomachinery, which tends to reduce E tr to a low 

value. According to Gruscbwitz (reference 6), observed bou n d ar y 
layers were always laminar if E < 250 regardless of the intensity 
of free-stream turbulence, and turbulent boundary layers have been 
observed for E > 250. Gruschwitz neglects to mention the pressure 
gradient at which these observations were made and as a consequence 
the assumption E^. r = 250 might be somewhat low for flows with zero 

pressure gradient. It should be noted that if for any reason transi 
tion is delayed so that the rapid diffusion rate prescribed for a 
turbulent boundary layer is applied to a laminar layer, separation 
results. On the other hand, the effect of pressure gradient is to 
reduce markedly the value of R-(- r , so that even should the transi- 
tion be delayed in the zero pres sure- gradient region it will prob- 
ably occur at the inception of the diffusion process. 


Turbulent Boundary Layer 
The momentum equation for the boundary layer is 


d§ 

dz 


+ u 2 (2+H) - la 


where 

t surface shear stress 

p gas density 

If the definitions for friction coefficient 




dU 2 


(9) 


(10) 


generalized momentum thickness 



( 11 ) 


and diffusion coefficient 


r. = I du 
a U dx 


( 12 ) 
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are substituted into the momentum equation as in reference 3, there 
results 


ar _ 7 
dx ~ 6 



(13) 


In order to complete the solutions of this equation for arbitrary 
Dp, tvo more equations are required: one to determine values of f , 

and one to determine the development of the form factor H. 


Region of constant velocity. - In the region of constant veloc- 
ity Dp = 0, so that the momentum equation reduces to 


dT 7 
dx “ 6 1 


(14) 


Falkner ' s data (reference 7) for surface shear vith constant veloc- 
ity show that 


f = 0.006534 (15) 

Equation (14) may therefore be integrated to give 

T - T tr = | f (x - Hr ) (16) 

where the subscript tr indicates the value at transition. 

Conditions for diffusion without separation. - The system of 
equations proposed by Gamer (reference 3) was selected to describe 
the development of the turbulent boundary layer because of its 
simplicity. In selecting a law for friction, Garner examines data 
for R > 1250 and shows that f is constant with a value of 
0.006534 if -0.01 < Dp < 0.004. (This is Falkner* s value of f 

for Dp = 0, reference 7). For Dp < -0.01 no correlation is 

obtained. In consideration of the lack of correlated data for high 
diffusion rates. Dp > -0.01 was selected for the range of dif- 
fusion rate and therefore f = 0.006534 is used for the friction 
coefficient. 

The development of the form factor H required for solution 
of the momentum equation is described by an empirical equation 
developed by Garner in reference 3. The relation is 
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T g = -e 5 ^ H " 1 ’ 4 ^ jDj, + A (H - 1.4)] (17) 

where A = 0.0135. 

The problem is then to define a velocity distribution that 
will avoid separation of the turbulent boundary layer. Experience 
shows that when H reaches a value of about 2.0, with conventional 
velocity distributions it subsequently rises rapidly and separation 
occurs. Values of H for separation are therefore usually given 
in the range 2.0 < H < 2.6. It is therefore proposed that one 
possible criterion for prescribing the velocity is that H rise 
to some value and thereafter remain constant and less than 2.6. 

Equation (17) describing the development of E shows that if, 
H < 1.4 + (-Bp/A) , then dH/dx > 0 and H tends toward 

1.4 + (-Dp/A). If H >1.4 + (-Bp/A), then equation (17) shows 

that dH/dx < 0 and therefore H always tends toward the value 
1.4 + (-Bp/A). 

If constant values of Bp are considered as possible condi- 
tions for determining the velocity distribution, then the momentum 
equation shows that for a given velocity change Bp should be as 

large a negative number as possible in order to reduce the surface 
for skin friction. (The rapid diffusion results in airfoils with 
relatively longer regions of constant maximum velocity and there- 
fore larger airfoil circulations.) For reliability of the boundary- 
layer equations, it is desirable to maintain Dp > -0.01. Therefore 

a possible condition is Bp = -0.01, the largest rate of diffusion 

with f = 0.006534. The limiting value of H for this case is 
2.14, a value that is probably, but not certainly, safe. Two values 
were therefore used in the velocity computations: Dp = -0.01, and 

Bp - -0 .006, which corresponds to a limiting value of H = 1.844. 

A third condition for prescribing the diffusion rate was also 
employed. According to reference 4, separation will result if 


— . — 0 dU _-0.08 T dU -0.247 
^S“Udx R = U dx 


BpE 


-0.247 


(18) 
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attains the value of -0.0013. The third type of velocity distri- 
bution computed was for the maximum diffusion rate with both 
Dp > -0.01 and Dg > -0.00115, a value arbitrarily assigned vith 
the condition that it be greater than the separation value -0.0013. 


Velocity distributions in the region of diffusion. - For the 
cases where Dp isconstant, equations' (12) and (17) are used to 

eliminate the differential of x, and to obtain 


dU _ % -5(1 - 1.4) dH 

U “ " A H - 1.4 + %/A 

which integrates to 


(19) 


U % 5Dp/A 


106 if •- 


|li ( -51-5I>p/A + 7) - Ei(-5I>p/A)J (20) 


where 



is the exponential integral function (tabulated in reference 8), 
having the principal part of the integral for y > 0 (that is, for 
H < 1.4 - Dp/A). The constant of integration for U *= Uq was 
chosen at 1 = Hq = 1.4. This value of H, U 0 , and Tq will not 
be included on any point of the airfoil if laminar diffusion pre- 
cedes turbulent diffusion. 


Equations (12) and (13) give 


dT 

T 


7 

6 



dU 

U 


( 21 ) 


For the factor dU/U of the second term on the right side, a sub- 
stitution is made from equation (19) to obtain an equation that 
may be integrated to 



_f_ 

% 




, _7_ Jr 
+ 30 A 


" -5(1 - 1.4) 

1-0 


( 22 ) 
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The arc length is then given by 



x 


n u 



T_ dU 
T 0 U 


(23) 


The curves for the boundary- layer parameters In the turbulent- 
diffusion region are shown in figure 4 for Dj = -0.010, Dp = -0.006, 

and maximum diffusion rate with Dp > -0.01 and Dg > -0.00115. 

In solving the boundary- layer equations (13) and (17) (with 
f = 0.006534) with the condition Dg > -0.00115, the value of Dp 

rises with B. The limitation Dp > -0.01 was added at the high 
Beynolds numbers in order to keep Dp in the range for which the 

friction law f = 0.006534 is valid. The solution of the equations 
in the region of variable Dp was obtained by approximating each 

portion of the curve with the solution for Dp = constant = mean 
value of Dp in that step. 


Eolations between Velocity Distribution 
and Blade Beynolds Number 

The relations between the velocity distribution and blade 
Beynolds number U f X/v are developed. Some of the results are 

given in graphical form, which permits the use of the generalized 
velocity-distribution curves developed for the laminar and tur- 
bulent boundary layers (figs. 2 and 4). These curves will give 
the velocity distribution in the region of diffusion if the scale 
constants Lj, Uq, and Tq are knownj the remainder of the sur- 
face is determined by the constant Xj or x t . 

As assistance in finding the appropriate constants, the quan- 
tity Xj/x (or x t /X for completely turbulent diffusion) is 

plotted in figure 5 as a function of the Beynolds number UfX/v 

and the velocity ratio U_/U . (The method of computing fig. 5 

x max 

is presented in appendix B.) The curves should be applied in the 
following manner: 

The Beynolds number U^X/v and velocity ratio Up/t^ nax are 
selected. Examination of figure 5 will indicate whether the 
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diffusion is laminar, mixed, or turbulent. If the diffvision is 
laminar, figure 5 gives the value of Xj, and then Lj = 0.441 x^, 
vhich establishes the scale for the use of figure 2. 

If the diffusion is altogether turbulent, figure 5 determines 
Xj. and then 



where (X - Xq)/T q is given in figure 4 as a function of U f /U 0 
and U 0 = 1^. 

The case of mixed diffusion is considered to be of no technical 
significance in this application because of uncertainty in locating 
the position of transition. The effect of Reynolds number is there- 
fore considered only in regions of purely laminar or purely turbu- 
lent diffusion, although x^/X and 0 /UpX were computed for 

®tr = 250 and Dp = -0.01 in the mixed-diffusion region. 

Suction-surface potential. - A further step in computing air- 
foil performance is the integration of the velocity curves to obtain 



a dimensionless form of the potential on the suction surface. 

Results of the computations are shown in figure 6. The comparative 
magnitudes of the potentials for various values of Dp (or Dg), 

Uf/ Umax * Keynolds number are significant because in .the 

system used for estimating blade circulation, an increase in suction- 
surface potential is generally accompanied by an increase in blade 
circulation. 

The obtainable potential is highest with the greatest diffusion 
rate (compare figs. 6(a) and 6(b)) because the Tnavirmim velocity is 
maintained along the largest portion of the airfoil surface. The 
most pronounced effect of Reynolds number is the increase in poten- 
tial shown in the change from laminar diffusion to turbulent 
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diffusion (fig. 6(a)). The increase in potential begins at a 
Reynolds number in the vicinity of 80,000 to 90,000 and for 
Uf/Ujnax > 0.6 is fully realized at the blade Reynolds number 

approximately 260,000 for Dg > -0.00115 and Dp = -0.006 and 
approximately 180, 000 for Dp = -0.010. This change Increases with 
decreasing diffusion-velocity ratio Uf/ Umav Although there is 

no effect of Reynolds number on suction- surface potential for dif- 
fusion with a completely laminar boundary layer, the effect of 
increase in Reynolds number for turbulent diffusion is a decrease 
of potential for constant values of Dp. An insignificant effect is 
shown, however, when the criterion of Kalikhman for separation is 
applied (fig. 6(c)). Because at large Reynolds numbers the condi- 
tions % = -0.00115 is replaced by Dp = -0.01, the curves so 

constructed approach the curves for Dp = -0.010 in that region. 

(The large Reynolds number condition is equivalent to the assumption 
that the extent of the laminar boundary layer approaches zero . ) 

The case of infinite Reynolds number is of particular signi- 
ficance because it is the most conservative design with turbulent 
boundary layer for constant diffusion coefficient Dp. The poten- 
tial at large Reynolds numbers was therefore computed and is shown 
in figure 7 for values of Dp = -0.010 and -0.006. The curve for 
Dp = -0 .010 at high Reynolds numbers is practically the same as 
that of Dg > -0.00115 for all Reynolds numbers because the poten- 
tial for the Dg-condltion shows so little variation with Reynolds 
numbers. 

The trend of curves shows that highest potential occurs at a 
UfAW of approximately 0.58 with very little variation between 

0.53 and 0.63. Therefore, if a design is considered with a pre- 
assigned value for the velocity diagram and U f , there is negligible 
gain in using higher peak velocities than U mpy ~ 1.6 U f . For dif- 
fusion with a laminar boundary layer, figure 3 shows negligible 
increase in obtainable potential for peak velocities U higher 
than 1.25 U f . 

Deviations from design conditions. - A knowledge of the range 
of operation of any cascade of airfoils is essential in determining 
its suitability for application. For the sets of blades considered 
here, the effect of change of Reynolds number is nonexistent with 
purely laminar diffusion. A very slight and indeterminate effect 
is shown for turbulent diffusion with Dg > -0.00115. If the limi- 
tations imposed by constant Dp are closer to the physical facts, 

then less lift is obtainable at higher Reynolds numbers. 
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Consequently, increase in Reynolds number of operation would 
involve the danger of separation. A safe estimate of lift would 
therefore he made on the basis of the potential for the maximum 
Reynolds number of operation provided the minimum Reynolds number 
of operation does not extend into the regions requiring laminar 
diffusion rates. If the maximum Reynolds number of operation is 
taken to he three times the minimum, then a very slight effect of 
design 1 Reynolds number on potential is indicated even for the dif- 
fusion rates Dp = -0.010 and -0.006. The Reynolds numbers at 

which some diffusion with laminar boundary layer occurs depend on 
the value selected for R tr (in this case 250), on the diffusion 

rate, and on the velocity ratio U^y /Up . Estimates of lift made 
on the basis of Dp = -0.010 and RN-»w would be in agreement 
with estimates based on Ig > -0.00115. If a blade designed on 
such a basis was used at low enough Reynolds numbers, then the 
laminar boundary layer existing in the region of diffusion would 
separate because of too steep a pressure gradient. 

If a rough finish or a higher angle of attack than the design 
value occur in application of the blade, an Increase in the value 
of T t over the design value may be expected. Then the initial 

value of Dp = T Qj will increase in the same proportion at 

the beginning of diffusion. Sample computations show, however, 
that for velocity distributions based on constant Dp, this 

initial increased value for Dp is not maintained but drops to 

some value intermediate between the initial value and the value 
for which the velocity distribution was prescribed. If the value 
of the form factor for separation E is assumed to be 2.2, it 
would correspond to the limit for a constant value of Dr given 
by 


Dp = - 0.0135 (H - 1.4) = - 0.0108 

If the design velocity distribution corresponded to 
Dp = -0.006 (limiting H = 1.85), and then in operation T were 

increased by 80 percent to give an initial value of Dp = -0.0108, 
the resulting limiting value of Dp would be greater than -0.0108 

(less in absolute value) and the limiting value of H would be 
less than 2.2; thus, according to this criterion, separation would 
be avoided. If the design velocity distribution corresponded to 
Dp = -0.010, however, an increase of 80 percent of T would result 

in separation according to this criterion. 
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ESTIMATE OF BLADE CIRCULATION AND SPACING 
Effect of Thickness 

An approximate calculation is made for the blade circulation 
from the upper-surface velocity distribution by a modification of 
the method, for thin airfoils. The basic idea is that the airfoil- 
surface velocity is the sum of two components, one of vhich results 
from the basic thickness form uncambered, and the second of vhich 
results from the curvature of the camber line. The suction-surface 
velocity of the thin, slightly cambered, isolated airfoil is then 
given by 

U' = U» + AU' - V* |; + 

vhere 

U' suction- surface velocity of airfoil 

t 

U' velocity component resulting from uncambered thickness dis- 
tribution 

AU' velocity increment resulting from curvature of camber line 
V' free- stream velocity 

The ratios U'/V' and AU'/V' are independent of the magnitude 
of V* . A similar equation holds for the pressure surface. 



U p ’ = U' - AU' = 2U' - U’ 


vhere U p ’ is the velocity on the pressure surface. If the 

suction- surface velocity is knovn, and if the velocity component U' 
resulting from the thickness can be determined, then the pressure- 
surface velocity can be computed. 


In applying these ideas to the cascade, the principal change 
is the use of the distorted flov field in vhich the airfoil is 
immersed; the quantity V* is nov regarded as a variable quantity. 
Near the leading edge, the velocity is substituted for V’ 

giving 




( 26 ) 
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where_the subscript 1 Indicates the value near the leading edge, 
and. U^'/v* is the velocity rise over the free-stream value for 

the isolated airfoil. The value for U’/V* is a function of the 
assumed thickness distribution. Similarly, near the trailing edge 
the free-stream value vised is V 2 giving 



(27) 


At the point of the airfoil where the air has been turned 
halfway so that the flow direction is given by the angle where 

tan ctjn = i (tan + tan aj>) (28) 


where 

a angle between velocity vector and normal to cascade axis 
Subscripts : 

1 value for upstream of cascade 

2 value for downstream of cascade 

the same procedure might be used as for the leading and trailing 
sections of the airfoil with the substitution of V_ for V' in the 

__ HI 

equation for U m . The value of V m is computed from the continuity 
equation 

cos o Vg cos a*, * V m cos (29) 

which does not include the effect of airfoil thickness in blocking 
the flow area and increasing the flow velocity. If the spacing s 
is not too small, the correction for the distorted flow that is 
applied is 



This formula is derived in appendix C. Then, 
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When the blade spacing is very small, however, It becomes more 
accurate to regard the flow between the blades as flow in a channel. 
From this standpoint, the flow velocity is increased over V m by 

the ratio of the flow area without blade thickness to the flow area 
with blade thickness. That is, the effect of thickness is given by 

U = V ( - 'N (31) 

• ■ Hf - v sec °v 

where t m is the blade thickness at the point x^ where the air 
flow has been turned to the direction a^. This equation is applied 

only if the point of application of the lift force x^ is within a 
range of values specified by 

| < x m < X - | (32) 


This condition was selected because it is believed that inside this 
range the effect of variations of inflow and outflow direction on 
the velocity is expected to be small enough to be neglected for the 
purpose of the calculation. Thus, the flow is expected to be 
essentially that inside a long channel. 


Relations have been developed for computing U at three points 
on the airfoil. The family of airfoil cascades being considered 
will be restricted by assuming a second-degree parabola connecting 
the points (0, U^), (x^ U m ), and (X, TJ f ). That is, 


- 2 - = — + (- 
uj uj ® X + B [l. 


(33) 


where 




U* 


u* 


„(2) U i/ U f ^m/ U f u f/ u f 

B = - (x m /x)d + Ti""^/xy 


Uf/Uf 


3' 


} ( 34 ) 
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Estimate of Circulation and Spacing 
x __ * _ 

By assuming that U*, and Uf are known, the curve for 

U can then he determined when the value of is found. Because 

is the center of pressure, it may he found as the weighted mean 
value of x; the weight assigned to each value is the local hlade 
loading. 


Because 


*m 


is at the center of pressure. 



(35) 


where 

cp potential on upper surface measured from leading edge to a 
point x 

r airfoil circulation, ($ - i> p ) 

$ value of cp at trailing edge 
The quantity, 

may he found from 



either figures 3, 6, or 7. The quantity. 


<t> 



xdcp = / xU dx = UfX 2 m 


may he found from figure 8. The curves for Bp = - 0.010 are 
approximately correct for Dg > -0.00115. 

For calculation of the potential cpp, the high camber of 

airfoils in cascade must he considered because it results in a 
substantially different length of the suction and pressure 
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surfaces unlike isolated airfoils, which are nearly straight. The 
effect of the difference in length is approximated by setting 

Xp = px (36) 

where p is some constant nearly equal to unity. The circulation 
is thus. 


U dx 


/ U p dxp * / ju - P (2U - U)] 
JO JO 


r - * - 2p U dx - P4> = $ - 4>, 


where 


4> p = 2P U dx - p$ 


Using equation (33) for U, there are obtained. 


l/ U f U i\ b (2) 


0 * ■ V - t t ) - 


xd<p p = p I x U p dx * p / x (2U - U) dx 


with 


_2 /l U i B* 1 ) *(2)' 

*[zir + V + V 
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The factor 3 Is estimated on the assumption that the airfoil 
suction surface and pressure surface are circular arcs intersecting 
at the leading and trailing edges. The angle of turning X of the 
mean line is assumed equal to the air deflection: 

X = - ct 2 

The turning X of the suction surface is then computed from 


sin \r. Js = 1 X cos 3/4 
4 X 4 sin (X/4) 


(40a) 


Then the turning of the pressure surface X^ is computed from 

tan = tan t - 2 ^ ^ 


4 X 2 sin 

2 


(40b) 


The ratio of lengths is 


_ jp \> sin X/2 
P- X “ X sin Xp/2 


(41) 


Substituting the values for the integrals (38) and (39) in 
equation (35), there results 


C (3) (h'f * c (2) (Osf * >> u » ~ °i 

C VX/ +C \X/ C V X/ 12 XL 


(42) 


where 


c (3) „ U 1 * °f . J_ 




2 U* 


2 P ) U f X 


c (2) _ /i±j\ j_ + . i % . s ^ 


2P J UjX 


2P 


3 Uf 6 U f 


,(1) = 5_ ^ 1 . 1+J. 

/ — m T J rr A O “4 TT 


r 12 U- 4 Uf. 2 p 


6 U* 


J 


( 43 ) 
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With the values for Uj., t^, and U f equation (42) may he solved 
for equation (34) for b( 1) and b( 2 ), and finally the curve 

for U (equation (33)) determined. Furthermore, equations (37a) 
and (38) determine the potential $> p , so that the dimensionless 
circulation T/u^X is then 


r . $ $p 

U f X = U f X " U f X 

For potential flow. 


(44) 


r = s (V^ sin - Vg sin a 

The velocity may he eliminated'by the continuity equation (24) 

resulting in 


X Vg cos 02 (tan - tan a^) 

s “ C X 

UjX 


(45) 


which determines the hlade spacing. 


Comparison with Exact Procedure 


For the purpose of Judging the accuracy of the method of 
estimating the pressure- surf ace- velocity distribution, a comparison 
was made with velocities computed for potential flow hy an exact 
procedure. The data were obtained from reference 9. The comparison 
is shown in figures 9(a) and 9(h). In both cases, values of U’/V' 
were selected to fit approximately the NACA 65-series isolated air- 
foil data (reference 10), which gave 


u i' ^ 

V7-- 1-025 


V 


1 + 1.35 ~ 


V 


> (46) 


J 


0.950 
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By' using_these values_and equations (26), (27), (29), and (30), 
Ui/V m , Uf/v m , and were computed. Next p was obtained 

By use of equations (40) and (41). Evaluation of the data for 
suet ion- surf ace- velocity distribution yielded and m(Uf/V m ) . 

It was now possible to substitute all these values into equa- 
tions (43 ) and solve for %/X by use of equation (42 ) . The solu- 
tion was then continued in the manner indicated in the previous 
section. 

Pressure- surface velocity was then plotted using the equation 


The shapes of the estimated and exact velocity distribution 
on the pressure surface are quite similar. The solidity computed 
from the estimated velocities is, however, 6-percent higher for 
figure 9(a) and 4-percent lower for figure 9(b) than the exact 
value. These values then give an indication of the order of _ 
inaccuracy of the method of estimating solidity. The error in 
pres sure- surface potential was 6 percent for the airfoil of 

figure 9(a) and 2 percent for the airfoil of figure 9(b). 


EFFECT OF SEVERAL PARAMETERS ON REQUIRED CASCADE SOLIDITY 

The relations developed in the previous section are applied 
to study the effect of several design and operating parameters on 
required cascade solidity. 


Computations 

The data given in the problems solved consisted of the fol- 
lowing: (1) Reynolds number U f X/ v (either very high or very 

low); (2) diffusion coefficient Dj. (either -0.010 or -0.006 for 
a turbulent boundary layer); (3) trailing-edge loading (U^-UfJ/Vg 

(assumed to be 0.1 because cascade data Indicated this value as 
attainable); (4) velocity ratio v /Vg ; (5) velocity- diagram 

components V^, Vg, o^, and ctg from which V m and are 

computed; and (6) velocity ratios U^’/V' , l^'/V', and Uf ' /V ' , 

which are determined from the assumed basic thickness distribution. 
In all cases the assumed values were U^'/V' = 1.025 and 
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U f /V‘ = 0.95, but U^/V' was varied to obtain various values for 

the thickness. All the assumed quantities were varied in the 
examples to show the effect on the required cascade solidity, the 
velocity distribution, and the thickness. The sequence of com- 
putation is: 

(1) Compute 



1.025 



0.95 




°f - U f \ 
v V 2 / 


1.05 



(2) From Uf/ U mAV , the assumption of either laminar or tur- 

bulent diffusion, and the desired rate of diffusion, figures 2 to 8 
are used to find the suction-surface-velocity distribution and the 
values for and m. 

(3) Values for the solidity X/s and the pressure center Xj„/x 

are assumed. The thickness is then computed from equations (30) 
and (46) if either %/s < 0.5 or %/s > X/s - 0.5. If 

0.5 < %/s < X/s - 0.5, then equation (31) is used. The factor p 
is then obtained from equations (40) and (41). 

(4) The quantities C< 3 4 5 ), C< 2) , C^ lj , %, B^ 1 ), B* 2) , ®/U 
r/U f X, and X/s are next confuted from equations (43), (42), (34), 
(38), (44), and (45). 

(5) Steps (3) and (4) must be repeated if the solidity X/s 
varies considerably from the assumed value or if, because of an 
incorrect assumption for %, an incorrect choice was made in the 
equation for the computation of %/x and (3. 
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Besults 

The results of the computations are shown in the following 
table and in figure 10. The table gives the assumed values for the 
design parameters and the computed solidity, whereas figure 10 shows ‘ 
the details of the velocity distributions with sketches to indicate 
the cascade solidity and the airfoil thickness. The circles drawn 
on the camber lines indicate the location of Xq and the airfoil 
thickness tjj. 


Case 

Boundary 

layer 

Diffusion 

coefficient 

*1 

a 2 

U max 

U f -U f 

*m 

Z 


(deg) 

(deg) 

• v 2 

V 2 

z 

8 

a 

Turbulent 

-0.010 

45 

38 

1.615 

0.10 

0.10 

0.30 

b 

Turbulent 

- 

.010 

45 

38 

1.312 

.10 

.10 

.64 

c 

Turbulent 

- 

.010 

45 

38 

1.312 

.10 

.06 

.50 

d 

Turbulent 

- 

.010 

45 

38 

1.312 

.00 

.06 

.65 

e 

Turbulent 

- 

.010 

45 

-45 

1.615 

.10 

.20 

1.6 

f 

Turbulent 

- 

.010 

45 

-45 

1.615 

.10 

.15 

1.3 

6 

Turbulent 

- 

.010 

45 

-45 

1.312 

.10 

.11 

1.8 

h 

Turbulent 

- 

.010 

0 

-45 

1.312 

.10 

.18 

1.0 

' i 

Turbulent 

- 

.006 

45 

38 

1.615 

.10 

.10 

.42 

j 

Turbulent 

- 

.006 

45 

38 

1.312 

.10 

.06 

.56 

k 

rm-r 

- 

.06618 

45 

38 

1.615 

.10 

.06 

1.6 

1 

T^nrt Tmr 

- 

.06618 

45 

-45 

1.312 

.10 

.08 

3.1 

m 

T^nvlmr 

- 

.06618 

0 

-45 

1.312 

.10 

.09 

1.7 


The effect of variation of maximum surface velocity on 

solidity Z/s of a compressor- type blade is shown by comparison of 
cases (a) and (b) where all design parameters are the same except 
VtAp* The blade with a lower value of U maT / V 2 is obviously 

better suited for operation with a compressible fluid with a high 
Mach number, but has the disadvantage of requiring a higher solidity 
(twice as many blades). 

If in case (b) the thickness is reduced from 0.10 to 0.06, the 
cascade corresponding to case (c) results, having only 78 percent of 
the solidity of case (b). This decrease results from the higher 
aerodynamic loading of the thinner blade because of the lower pressure 
surface velocities caused by the thickness. 

In cases (a), (b), and (c), the trailing- edge loading corre- 
sponds to an air flow that does not satisfy the Kutta condition. 

The actual blade having the prescribed velocity distribution would 
have an extension from the point x/Z = 1.0 to larger values in 
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which the loading would drop to zero. Similarly, the initial parts 
of the velocity distributions do not include the leading edge. 

That is, the velocity distributions shown do not extend over the 
entire airfoil, but do include practically all the aerodynamically 
loaded region. A velocity distribution was computed for zero trailing- 
edge loading (case (d)) in order to observe what sacrifice would be 
involved in avoiding separation at the trailing edge and attempting 
an airfoil design that would satisfy the Kutta condition. Compari- 
son with case (c) shows that the solidity is about 1.3 times as 
great. The effect would be less in the case of reaction blading 
where the inlet velocity results in large loading in the initial 
part of the airfoil and a relatively smaller effect of trailing- 
edge loading. The effect of decrease in t railing-edge loading is 
also less for blades of small diffusion ratio lI maT /U-f>, because in 

such cases the suction-surface velocity is maintained at its maximum 
value to a point very close to the trailing edge and, as a consequence, 
the loading is also maintained over most of the airfoil. 

For an impulse blade (V 2 = V^) of large turning (turbine type), 

the velocity distribution is very much the same as for a compressor 
blade with nearly the same velocity curves U and U, although the 
thickness will radically vary because of the low value for V m . For 
example, compare cases (a) and (e). 

If the turbine-blade thickness is reduced to the value closer 
to the compressor blade, the circulation is increased (case (f ).) , as 
might be expected. The large diffusion ratio on the pressure surface, 
which practically insures the existence of a separation bubble there, 
should also be noted. 

The effect of decreasing the maximum surface velocity % is 

shown in case (g), which has lower blade circulation even though the 
thickness was decreased from 0.15 to 0.11. Solidity increased 
39 percent over the value for case (f). 

For an inlet guide vane turning the air from an angle of 0° to 
an angle of -45° (case (h) ) , the velocity distribution changes from 
the impulse cascade (case (g)) because of the large decrease in the 
inlet velocity. Sven with the increased thickness, local velocities 
on the pressure surface do not involve high diffusion rates because 
of the over-all Increase in velocity (V^/Vg = 0.707). 

The effect of reducing the rate of diffusion from that of 
cases (a), and (c) (lip = -0.010) to a lower rate (Up = -0.006), 
while maintaining the same values of the other design parameters, was 
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determined by computing the corresponding cases (i) and (j), respec- 
tively. The positions at which the U/V 2 values axe equal for the 
two cases are where the values of a are the same. There is some 
chordwise shift in the location of a given value of a when changing 
from case (a) to case (i) because of the shift in the loading dis- 
tribution. Consequently, the location of shifts and the thick- 

ness distribution changes. The last half of the blade (case (i)) 
does very little turning so that the solidity increases from 0.31 to 
0.42. The usefulness of this part of the blade could be greatly 
increased if a larger turning of the air were involved, for the 
U/V 2 curve would rise to higher values in the first half of the 

airfoil and drop to lower values in the last half, thus shifting the 
load backwards on the blade surface. The change in velocity and 
solidity from case (c) to case ( j) is slight because the rate of 
diffusion has very little effect on suction-surface potential when 
the diffusion is small. (See fig. 7.) 

An effect of blade Reynolds number UfX/v is shown by compari- 
son of cases (a) and (k). In case (k), the expansion is entirely 
laminar (low Reynolds number), whereas in case (a) the extent of the 
laminar region is assumed zero (high Reynolds number) . Even though 
the laminar-layer blade is only 6-percent thick, the solidity (1.6) 
required is five times as large as for the blade with a turbulent 
boundary layer. This loss in circulation results from the fact that 
laminar diffusion must be begun so much earlier when the boundary 
layer is thinner. On the pressure surface of the blade is to be 
observed a region near the trailing edge where the rate of diffusion 
increases to a value that will certainly induce separation. A better 
aerodynamic design could be obtained with a thinner blade, which 
would be necessary if the turning were greater. Some improvement 
would also probably result with, a less simple mean velocity curve. 

It appears safe to conclude that it is extremely difficult, if not 
impossible, to obtain a compressor blade of reasonable thickness and 
turning, if flow separation is to be avoided with a laminar boundary 
layer. 

The design of an impulse blade with laminar boundary layer 
appears to be entirely feasible (case (Z)) as a consequence of the 
low value for U m /v 2 , which is a result of the low value for 

Vjjj/Vr. A reaction blade with laminar boundary layer also is fea- 
sible. (See case (m).) The percentage increase in solidity for 
cases (Z) and (m) over that for cases (g) amd (h) is much less than 
the percentage increase in solidity for case (k) over that for 
case (a). As a general rule, cascade solidities for laminar and 
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turbulent diffusion can be expected to be more nearly alike vhen 
diffusion is small (Pma-r/Uf ~»1.0) than when diffusion is large 

(W u f^ 2 -°)> because of the smaller difference in suction- 
surface potential. (See fig. 6(a).) 

COMPARISON OF EXPERIMENTAL WITH INDICATED ATTAINABLE 
' SOLIDITY VALUES 

In order to estimate the degree with which blades in use approach 
the maximum attainable circulation as indicated by the present method 
of analysis, cascade solidities were computed for conditions com- 
parable to those for which performance data were available. It was 
desired to compare the solidities of the tested blades and those 
corresponding to the proposed velocity distributions on the basis of 
operation at very low Mach numbers, very high Reynolds number, the 
same values for a^, cl,, maximum thickness t/X, Umav/ Vp^ and 

Uf/V 2 . The data available did not give all this information, or did 

so with no indication as to the efficiency of the cascade. Data used 
were obtained from references 11 and 12. In none of the comparisons 
shown is it possible to evaluate cascade efficiency, although for 
every case the comparison was made on the basis of the optimum oper- 
ating condition of the cascade as stated. The comparisons are shown in 
the following table: 


Refer- 

ence 

Blade 

desig- 

nation 

“1 

(deg) 

“2 

(deg) 

^max/^2 

Exper- 

imen- 

tal 

X/s 

Estimated x/s , 

Dp = -0.010 

Dp = -0.006 

ii 

65- (18)10 

45 

15 

1.699 

1.063 

1.68 a 


11 

65- (12)10 

45 

23.8 

1.602 

1.034 

1.17 

1.60 

11 

65- (4)10 

45 

37.0 



.78 

.88 

11 

65- (0)10 

45 

42.9 

1.186 

1.003 

.36 

.38 

12 

64-(A)06 

0 

-45.8 

1.180 

1.538 

1.12 


12 

64-(B)06 

0 

-52.0 

1.133 

1.534 

1.20 



a t/x = 0.0664. 


All of the cases are directly comparable as previously described 
except the 65- (18)10 blade, for which it was impossible to find a 
solution with 10-percent thickness, indicating a strong probability 
of separated flow on the blade. A solidity was computed with the 
thickness reduced to 6.64 percent. 
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The 65- (12) 10 blade has a lower solidity and therefore a higher 
circulation than either of the estimates for a maximum- circulation 
blade. Possible causes for this discrepancy are: (l) slight flow 

separation from the surface of the 65^( 12 ) 10 blade, and (2) the 
assumption of a parabolic curve for U, which may involve some sac- 
rifice in attainable circulation. For lower diffusion (65- (4)10 and 
65-(0)10) and for reaction blades (64-(A)06 and 64-(B)06), the estimates 
indicate that higher circulations are attainable. The 64-(A)06 and 

64- (B)06 blades have such high values for U^/ U^t that there is no 

substantial difference when the lower diffusion rate Dp ■ - 0.006 

is used in the solidity estimate. It seems reasonable to expect that 
higher blade circulations are obtainable with blades other than the 

65- and 64-series for low pressure rise and for pressure drop. The 
65- (0)10 blade is particularly subject to improvement. 


CONCLUSIONS 

The analysis of the limitations on the circulation about blades 
in cascade indicates that under the assumption of the several criteria 
for separation and the avoidance of local separation of the flow, the 
following conclusions can be drawn: 

1. For a preassigned maximum velocity and suction-surface 
trailing-edge velocity, the suction-surface potential and airfoil 
circulation increases with the Increase in permissible diffusion rate. 

2. For a laminar boundary layer, Loitsianskii's equation deter- 
mines blade circulations independent of the Reynolds number. For 

a turbulent boundary layer, Kalikhman's criterion for safe diffusion 
(Dg > - 0.00115) indicates very slight changes in blade circulation 

with variation in Reynolds number. If a constant diffusion rate 

(Dp = - 0.010 or - 0.006) is used, however, some slight decrease 

✓ 

in attainable potential with increase in Reynolds number is indicated. 
There is a very marked increase in obtainable lifts when the design 
is changed from that for a laminar boundary layer to that for a 
turbulent boundary layer. This effect decreases with decrease in the 
diffusion ratio U^y/ Up or Um^ /Vg. With transition at momentum- 

thickness Reynolds number R = 250, the change begins to occur at a 
blade Reynolds number RN between 80,000 and 90,000. For completely 
turbulent diffusion with Uf/ Up^y >0.6, the blade Reynolds numbers 

are 260,000 (Kalikhman's criterion and D T = - 0.006) and 
180,000 (Dp = - 0.010). 

3. Suet ion- surface potential and consequently the airfoil cir- 
culation increases with increasing maximum velocity ratio or 
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U ma-r /U-p? however, negligible gain in circulation is obtained with 
Umbt /U-p greater than 1.6 with turbulent diffusion and greater than 
1.25 with laminar diffusion. 

4. An airfoil designed for laminar boundary layer may be operated 
at any Reynolds number. An airfoil designed for a turbulent boundary 
layer has a lower Reynolds number limit for operation without sepa- 
ration. This limit is a function of the velocity ratio Umav /U-p, 

the diffusion rate, and the condition for transition. For R = 250, 
and diffusion velocity ratios U m n T / Uf < 1*67, the lower limit of RN 

is equal to or less than 260,000 for small diffusion rates 
(Dp = -0.006, Dg = -0.00115) and equal to or less than 180,000 

for larger diffusion rates (&p = -0.010). 

5. Improvement in airfoil circulation always results from adding 
loading at the trailing edge, but the Improvement is less with 
decrease in over-all diffusion (decrease in V^/Vg) and with decrease 
in surface diffusion (decrease in U ^a-rA Jf ) ♦ 

6. Change in suction-surface-velocity distribution resulting 
from change in diffusion coefficient or character of the boundary 
layer has a large effect on blade circulation except when the dif- 
fusion is small. This small effect is to be expected in blades 
designed for high-speed operation and for reaction blades. 

7. For a fixed suction-surface-velocity distribution, obtainable 
blade circulation decreases with increasing blade thickness. 

8. The large discrepancy between the circulation of the 

65- (18)10 blade with the recommended velocity diagram and the esti- 
mated attainable circulation implies a strong possibility of sepa- 
ration of the flow from the 65 -(18) 10 blade. For the recommended 
conditions of operation of the 65- (4)10, 65- (0)10, 64-(A)06, and 
64-(B)06 blades, estimates indicate that higher blade circulations 
are attainable. 


Lewis Flight Propulsion Laboratory, 

Rational Advisory Committee for Aeronautics, 
Cleveland, Ohio, June 6, 1949. 
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APPENDIX A 


SYMBOLS 

The following symbols are used in this report: 


A 

0.0135 

B^,B^ 

coefficients of (x/x) and (x/x) 2 terms in 
r equation for U 

c (D,c(2), c (3) 

coefficients of (r^/X), (^/X ) 2 , and (x^X) 3 

terms in equation for ^ 

C 

chord 


diffusion coefficient with laminar boundary layer, 
/L dU\ 

\U dx/ 

B S 

(is®"**' 08 ) (according to Kalikhman, separation 
occurs when Do * - 0.0013) 


% 

Ei(y) 

F 

f 

H 

L 

m 


diffusion coefficient with turbulent boundary layer, 
(T dU\ 

\U dx/ 

exponential Integral, (I dk 

V-» k J 

parameter in approximate equations for laminar 
boundary layer 

coefficient of friction, R 1 / 6 ^ 

form factor, (6*/e) 

generalized momentum thickness of laminar boundary 
layer, (6R) 




xU dx 

o 

momentum- thickness Reynolds number, (U0/u) 
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blade Reynolds number, (U^X/u) 

spacing of airfoils along cascade axis 

generalized momentum thickness of turbulent boundary 
layer, (0R 1 / 6 ) 

thickness of airfoil 

suet Ion- surface velocity Just outside boundary layer 

arithmetic average of suction- and pres sure- surface 
velocities 

stream velocity 

V m corrected for thickness effect of airfoils In cascade 

total suet ion- surface length 

arc length on airfoil suction surface 

angle between velocity vector and normal to cascade 
axis 

ratio of pressure-surface length to suet ion- surf ace 
length, (ip/x) 

circulation around one airfoil 

displacement thickness of boundary layer 

✓ 

momentum thickness of boundary layer 
turning angle of upper surface 
turning angle of camber line 
kinematic viscosity of gas 
gas density 
surface shear stress 


velocity potential 
edge 


•a; 




complete velocity potential, [ f Udx 

Udx) measured from leading 



NACA TN 1941 


33 


Subscripts: 

f 

i 

1 
m 

max 

P 

t 

tr 

0 

1 

2 

Superscript: 


final value near trailing edge 

Initial value near leading edge 

value at beginning of laminar diffusion process 

value corresponding to mean of upstream and down- 
stream velocity vectors 

maximum value 

k* 

value on pressure surface 

value at beginning of turbulent diffusion process 
value at transition 

value corresponding to H = 1.400 in turbulent 
boundary layer 

upstream value 

downstream value 


pertaining to isolated airfoil 



34 


NACA TN 1941 


APPENDIX B 


METHOD OF COMPUTING FIGURE 5 

For diffusion with a laminar boundary layer, equations (7a) 
and (8) determine for any given velocity ratio Uf/U the value 
for xj. Then 

. Lj = 0.441 Xj 

For the case of diffusion with partly laminar and partly tur 
bulent boundary layer, a number of values vere selected for 


_ e l 

l 


the Reynolds number at the beginning of laminar diffusion. The 
units of velocity and length may be taken as Uj and Xj, 
respectively. Also 


Lj = 0.441 x z 


which establishes the scale of the laminar-diffusion region. At 


the transition point, Rj- r = 
is obtained from equation (6d) 


0 tr u tr 


u 


= L l l ”e7 


= 250, so that for L fcr there 

ZCF+Dl) 

R tr \F +ZD L 


and 


U H 


'tr ^tr 


U, 


max 


.R, 


2Dl 
F + ZI^ 


Because, 


L = 0R 
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then 


T tr “ e tr (®tr) 


The Initial condition for E la 


- , 4 . ^tr_ 

nr - 1#4 0.0135 


(reference 3), where 


%,tr= (§ g) tr ■ A<Ktr) 5/6 


For this value of figure 4 identifies U^/Uq and Rfc r /Ro> 

thus establishing the scale for the velocity U and the basic 
length Tq. By proceeding to the final velocity Up, the length 
is established, and KN = UpX/u may be computed. 

If the diffusion is entirely turbulent, then by assuming a value 
for the Reynolds number at the beginning of turbulent diffusion R^, 
there are obtained (xfc r and U^r are taken as units of length 
and velocity, respectively) 

Ar'0-^5* 


T tr ■ 


0.441 ^ 

5/6 

(Bte) 


Dj,,t r = 0 

Htr = 1.4 
U 0 = U max 
T 0 = T tr 
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APEEUDIX C 

EFFECT OF AIRFOIL THICKNESS ON FLOW THROUGH CASCADES 

By assuming that the effect of the airfoil circulation on the 
velocity at the point Xq, (where the flow has turned to the direc- 
tion otjg) is taken into account by the turning effect, the change 
in velocity from V]_ to Vj,, is given by 

V m = V x cos sec % 

The effect of thickness is now separately evaluated as the approxi- 
mate effect of a cascade of uncambered airfoils of stagger 
with no circulation. The thickness distribution is assumed similar 
to that of a Joukowsky airfoil, which is discussed generally in 
reference 13. The potential function W for a unif orm stream of 
velocity V flowing about a Joukowsky airfoil of small thickness 
can be shown to be 


W = V 


t * 


(1+2c)c 2 /16 

fT . 


where 

e real constant determining airfoil thickness 

. a2 

i complex parameter defined by z = £ + 


and the condition that 


with 



£ be a continuous, single- valued function 
where z is the complex position coordinate 


in the airfoil plane. The airfoil profile is given by the additional 

equation for a circle of radius (l+e)£ and center - cc/4: 

4 


t + f = ! (1+£ > el “ 


When the central angle of the circle to has a value of 2/3 n, the 
airfoil has a maximum thickness and at this point the velocity is 

U = V(l+2e) 
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If the potential function is expanded in powers of l/ £, the indi- 
vidual terns may then in turn be expanded in powers of l/z, giving 
the resultant approximation for large z 


w = y 



ec 2 /8 \ 
z+c/8 J 


This expansion neglects terns in l/z 4 and higher powers. From 
this expansion, the velocity at large distances from the airfoil is 


dW _ v _ Yc c 2 /8 _ _ (U-V)c 2 /16 

dz (z + c/8) 2 (z+c/8) 2 


At a half chord above the airfoil 
disturbance velocity — - 

16(z+c/8r 



the error in the 


is 26 percent; whereas the error 


in the entire velocity is (for U = 1.2 V) 0.8 percent. At a 
distance of 1 chord, the errors are 10 percent and 0.1 percent, 


respectively. At 1^ chord, the values are .5 percent and 0.02 percent. 


This approximate velocity distribution is seen from the f on of 
the velocity equation to be that of a doublet of strength 

2 

M = 2n(U-y>Ye 


placed in a uniform stream of velocity V. 

If a series of such doublets are placed on an axis with direc- 
tion ie” iam and with locations 

z 0 ± inse -10 ^ 11 ' (n = 0, 1, 2, . . . ) 
a potential function W will result with the equation 


w _M_ 
W “ 2* 


z - z 


+ Y" 


° \z - z Q + iru3e" ictm z - z 0 - inse" ic ^/ 
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or 


W = 


Me lq m 

21s 




If the potential due to the central doublet at z 0 is removed to 
find the modification of the uniform flow in which the central blade 
is located, there results 


W = 


M 


2ise 


-ion 


cot 


- &q) . 


ise 

ji(z 



By expanding in powers of (z - z Q ) and by neglecting powers higher 
than 2, there is obtained the approximation 


H ■ £ f B 2 e21%1 < 2 ' Z °> 

with the complex velocity 

dW . _M_ j£_ 2iajj, 

dz " 2« _ 2 
38 

When added to the mean flow, the component normal to the mean flow 
may be neglected, giving 

V ■ \ + & 75 008 2 °>n 

’ 3s 

The doublet strength has been shown to be 


M = (D - T m>0 ) = ^ (D - T m;9 ) 

If this equation is substituted in the equation for V m c , and the 
equation 

U U' 


used to eliminate U, there results 
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v m,c 


1 - fl ^ (£ - x ) 003 *“■ 


Or, if V_ „ ' is eliminated, there is obtained 

' ID j U 


u* 


= ^m,C y7 = 


U' l V' 
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Figure 2. - Development of velocity and laminar boundary layer with constant diffusion coefficient. 
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Figure 3. - Relation between suction- 

surface potential and velocity ratio 
for laminar diffusion with 0 L ■ -0.06618. 





Figure 4. - Turbulent boundary-layer growth. 





igure 4 . - Continued. Turbulent bounda ry-l ayer growth. 
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Figure 4. * Concluded. Turbulent boundary-layer growth. 
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X/*x pue x/^* 'sojjej qj6ua"| 


Figure 5. - Variation of length ratios x^/X and x t /X with blade Reynolds number. R tr * 250. 
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Blade Reynolds number, UfX/v 
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X 


igure 6. - Continued. Variation of suction-surface potential with blade Reynolds number. R tr ■ 250. 




igure 6. - Concluded. Variation of sue t i on-^su rface potential with blade Reynolds number. R tr * 250. 
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Figure 7. - Relation between suction- 
surface potential and velocity ratio 
for blade Reynolds number approaching 
infinity. 
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